On the Topological Origin of Entanglement in Ising Spin Glasses 
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The origin of thermal and quantum entanglement in a class of three-dimensional spin models, at 
low momenta, is traced to purely topological reasons. The establishment of the result is facilitated by 
the gauge principle which, when used in conjunction with the duality mapping of the spin models, 
enables us to recast them as lattice Chern-Simons gauge theories. The thermal and quantum 
entanglement measures are expressed in terms of the expectation values of Wilson lines, loops, and 
their generalisations. For continuous spins, these are known to yield the topological invariants of 
knots and links. For Ising-like models, they are expressible in terms of the topological invariants of 
three-manifolds obtained from finite group cohomology - the so-called Dijkgraaf-Witten invariants. 
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The history of modern physics is replete with stories 
of the progress that followed every time the notion of 
instantaneous action at a distance between two objects 
was sacrificed in favour of the concept of an interaction 
between the participants. With an appropriate choice of 
gauge fields acting as mediators of the interaction, we 
can account for most of what is known about the funda- 
mental forces of nature. Several examples readily come 
to mind viz. Maxwell's theory of electromagnetism, the 
standard model of the electroweak forces, quantum chro- 
modynamics - the theory of strong, nuclear forces - and 
Einstein's general theory of relativity for gravitational 
forces jjj. Drawing inspiration from these examples, we 
shall use the idea of gauge invariance as a tour de force 
to gain insight into the origin of entanglement in Ising 
spin systems. 

The need to understand the origin of quantum entan- 
glement arises from the recognition that it is the theo- 
retical bedrock that supports the emerging revolution in 
storing, processing, and retrieving information Q. Most 
of the work on the physics of quantum entanglement is 
- in view of its aforementioned technological import - 
focussed on finding ways to quantify it, in developing so- 
phisticated algorithms to handle computational complex- 
ity, in grappling with the issue of decoherence, and last 
but not least, in producing experimental devices which 
shape the promised revolution. Under the circumstances, 
very little time is invested in finding an answer to a fun- 
damental question: What is the origin of entanglement 
in quantum systems? We will show that the gauge prin- 
ciple once again holds the key to unlock this, yet another, 
mystery of nature. 

In view of the fact that entanglement is a purely quan- 
tum mechanical feature, unlike the fundamental forces 
of nature alluded to earlier, it may seem that the gauge 
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principle is inapplicable in this context. Such apprehen- 
sions are unfounded because it is by now well-known that 
purely quantum features with no classical analogues are 
often associated with non-trivial topological aspects of 
the theory. The simplest example which illustrates this 
intimate relation is that of a free nonrelativistic point 
particle moving on a circle. A total derivative term 
added to the Lagrangian of this model leaves the clas- 
sical mechanics of the particle unchanged, but leads to 
different, inequivalent quantizations of the particle due 
to the topology of the circle Q . Another example con- 
cerns identical particles which acquire anyonic quantum 
statistics in a plane, due to the non-trivial topology of 
the configuration space |4j. The important role played 
by topological objects like monopoles and vortices in pro- 
ducing confinement, and instantons in producing anoma- 
lies, in Yang-Mills gauge theories, is also well documented 
The reasons for the quantization of charge in elec- 
trodynamics, flux in superconductivity and conductance 
in the Hall effect can all be traced to their topological 
antecedents 6] . It should be noted that for each of these 
examples, there is no classical counterpart; the physical 
effects have a topological origin, and are most elegantly 
captured by the gauge principle. We then conjecture that 
quantum entanglement is also purely topological in origin 
and substantiate this claim by using the gauge principle, 
in the concrete example of the three-dimensional Ising 
model. 

We mention that an intriguing analogy between quan- 
tum entanglement and classical topology was first pub- 
licised by Arvind [3, following his serendipitous discov- 
ery of the similarity between the entanglement properties 
of the Greenberger-Horne-Zeilinger state || and the cu- 
rious linking properties of Borromean and Hopf rings. 
The analogy was further developed by Kauffman and 
Lomonaco 9] who examined the parallels between en- 
tanglement of states in quantum mechanics on the one 
hand, and the strings of a braid on the other. The hope 
that the breakthrough in finding new topological invari- 
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ants of knots and links |10| may be used to character- 
ize/classify entangled quantum states through the above 
analogy, gave rise to some initial excitement. The limita- 
tions of such a hope, as delineated in 0, were primarily 
because the much-needed rule to associate a closed loop 
to a quantum state, and the snipping of a loop to a mea- 
surement, could not be formalised. In the present work 
we would like to approach this problem from Heisen- 
berg's point of view, and deal directly with the physically 
observable quantities. This approach enables us to cir- 
cumvent the need to define the topological equivalents of 
spins and measurements. We thus overcome the tempta- 
tion to look for parallels between quantum entanglement 
and classical topology. Instead, we search for the intrinsic 
topological content of entanglement in a prototypical spin 
system namely, the three-dimensional Ising model. In 
the process, we discover unambiguous relations between 
the usual measures of entanglement 0] an d the expecta- 
tion values of gauge-invariant observables in a Z2 lattice 
Chcrn-Simons gauge theory. As is well-known, the lat- 
ter are just the observables used to calculate topological 
invariants of knots and links in topological field theories 
with continuous groups F° r the case of finite groups, 
the appropriate topological invariants, obtained by using 
finite grou p co homology, are called the Dijkgraaf-Witten 
invariants 

The Ising model is, as usual, defined by the Hamil- 
tonian H = — Jj2<ij> SiSj, where J > and Si = ±1. 
Here i, j label the sites of a three-dimensional cubical lat- 
tice and the < > parantheses indicate that the summation 
is over different, but nearest neighbour sites. The posi- 
tivity of J implies that the ferromagnetic state minimises 
the energy. For the opposite sign of J, the antiferromag- 
netic state would be favoured. As already motivated, we 
proceed to replace the above Hamiltonian which couples 
spins at different sites, by 

Hu = -jJ2 SiUijSj (1) 

<ij> 

in which the interaction between spatially separated spins 
is mediated by the gauge field U. The Uij live on the links 
connecting sites i and j, are ^-valued, and hence equal 
to ±1. The gauged Ising model presented above accom- 
modates either a ferromagnetic or an antiferromagnetic 
bond between various nearest neighbour sites 0. Such 
models play a crucial role in understanding the behaviour 
of disordered systems; in which subject they are referred 
to, and extensively studied, as spin glasses (l5|. In this 
letter, however, a spin glass is merely an expedient to re- 
place the action at a distance (of the order of the lattice 
spacing) between Ising spins by an interaction mediated 
by the gauge fields U^ . 

We define the partition function for the above model 
in the usual way as the trace over the Gibbs' measure 
i.e. Z = Tr e~ mu where /3 = l/k B T, k B being the 



Boltzmann constant, and T, the temperature. In the first 
step, called quenching in the language of spin glasses, we 
sum over all the spin degrees of freedom to get 

Z = A(J)J2^ Sv ~ Scs (2) 

u,v 

where 

S v = -JpJ2Ii v ii ( 3 ) 

□ □ 

and 

5 CS =/3^^(i-n^)(i-^) (4) 

<ij> □ 

The constants A and J are defined in terms of the 
strength of the magnetic interaction J, as follows: 

A( J) = coshJ and tanhJ = e~ 2J (5) 

The Vij , like the Uij , are ^-valued fields, but live on the 
links of the dual lattice obtained by a rigid translation of 
each site on the primary lattice by half a lattice spacing 
uniformly along each of the coordinate axes. The prod- 
uct of the dual gauge variables Vij around an elementary 
plaquette on the dual lattice is indicated by the □ un- 
der the product symbol. When the □ appears under the 
summation symbol, it is an instruction to sum over all 
such elementary plaquettes in the dual lattice. The Sy 
term is readily recognised as the standard Wilsonian ac- 
tion for the gauge field V. The Scs term is physically a 
measure of the flux passing through a dual plaquette per- 
pendicular to a given link U on the primary lattice. It is 
a Chern-Simons action on the lattice. Two crucial steps 
in arriving at the above result consist in an expansion of 
Z in the characters of the Zi group, and the introduction 
of the dual lattice variables. The details can be found in 
(lif. The generalisation to other finite abelian groups 
Z v follows along the same lines 0], as does the limit- 
ing case p — > 00. The latter corresponds to the familiar 
lattice U(l) Chern-Simons theory |18j . 

At this stage, a few general observations regarding the 
origin of the various terms in (2) are in order. First, from 
the continuum perspective, for the case of a continuous 
group, this result is easily anticipated. It is well-known 
from [19| that a derivative expansion of the fermionic de- 
terminant in this case, obtained by integrating out the 
fermions coupled to a background gauge field in three- 
dimensional spacetime, produces at the lowest two or- 
ders, the Chern-Simons and Maxwell terms. The result 
in (2-5) is a lattice realisation of the above celebrated 
continuum result. Second, from a purely lattice point of 
view, it is well-known that the three-dimensional Ising 
model is equivalent to a Z2 lattice gauge theory on the 
dual lattice j2jj. The Ising spins with nearest neighbour 
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interactions on the primary lattice have thus been traded 
for the ^-fields with a Wilson action that appears in (3). 
In the duality transformation, the U fields introduced 
by the gauge principle on the primary lattice are mere 
spectators and hence we arrive at a set of two ^-valued 
fields. Moreover, since there is an inextricable linkage 
between the primary and dual lattices, each link on the 
primary lattice pierces a plaquette of the dual lattice (and 
vice versa), and the Chern-Simons action is precisely a 
measure of this flux. Third, the noticeable difference be- 
tween the continuum and lattice realisations of what is 
essentially the same result is reminiscent of the lattice 
fermion doubling problem and has been discussed before 
in attempts to discretise Chern-Simons gauge theories 
[2l| . It may be recalled for the sake of completeness 
that, the mathematical reason for the doubling is that 
the Hodge star operator appearing in the Chern-Simons 
term couples cochains of a simplicial decomposition with 
the cochains of the dual decomposition. From a physi- 
cal point of view, the Chern-Simons term, owing to the 
fact that it has no independent gauge-invariant dynam- 
ics, couples matter fields to the magnetic flux. In the 
present case, the Ising spins residing on the sites of the 
primary lattice are the matter fields. The magnetic flux, 
in lattice gauge theory, is defined by the gauge fields on 
plaquettes. There is no natural way of coupling these two 
objects without doing violence to the structure of the lat- 
tice. This very fact was used by Kantor and Susskind for 
one of the early models for anyons |22^ . 

Finally, we mention that models of topologic al q uan- 
tum computation using anyons were proposed in |23j and 
experiments with flux-qubits |24j are beginning to realise 
some of these ideas. Quite independently, remarkable 
progress has been achieved by using statistical mechan- 
ical techniques to study spin glasses in theories of in- 
formation processing [23 following the proposal of |2fi| . 
Interestingly, equation (2) establishes a connection be- 
tween the above two seemingly different approaches to 
the subject of quantum information. More importantly, 
since the conceptual roots of the former approach lie in 
topology, while those of the latter lie in gauge invariance, 
it reinforces the idea anticipated at the beginning of this 
letter. 

To proceed further, we note that although equation (4) 
does not treat U and V on the same footing, the exponen- 
tial of —Scs which appears in equation (2) is invariant 
under an exchange of U and V |lfi| . However, equation 
(2) itself is lopsided because it does not have a Wilson 
term associated with the U field. This suggests that in 
the next step, called configurational averaging, where the 
quenched gauge degrees of freedom are summed over sub- 
ject to a given distribution, we choose the weights such 
that the U <-> V symmetry is restored. This is tanta- 
mount to using Hjj + Su where Sjj — ~KJ2a Yia Uij 
instead of Hjj. It is clear from the context that, in this 
case, the plaquettes under consideration belong to the 



primary lattice. If we now introduce the two-component 
vector Q = (U,V) and the matrices 

-(Si)- < 6 > 

the partition function can be rewritten in the neat form 

n 

where 

5 =§ E^(i-nw- fi )-^En^ ( g ) 

<ij> □ □ □ 

The < ij > in the above equation refers to, as before, 
nearest neighbour sites, but these could now be either on 
the primary or dual lattice. Next we take the infra-red 
(low momentum) limit and drop the TV-term to get, 

<ij> □ 

This is recognised as the lattice version of the topological 
Chern-Simons theory encountered in the studies of the 
double layer quantum Hall effect (2t| , and sometimes also 
referred to as the BF-theory |28j . 

We are now in a position to examine the entanglement 
properties of Ising spins in the above system. The ob- 
ject of central interest in studying entanglement is the 
(reduced) density matrix. Once it is known, a relatively 
straightforward calculation yields the von Neumann en- 
tropy S of the system through the standard formula 
S = — Trplnp. The single particle reduced density ma- 
trix pi, obtained by tracing over all the spins except the 
i-th spin, can be expanded as pi = | Y^ a =o c a (J i' where 
a — 1 and erf ^° are the Pauli matrices at site i. The re- 
quirement that Tr pi = 1, the reality of H , and the global 
spin-flip symmetry of the Hamiltonian imply that apart 
from Co which is unity, c\ is the only non-vanishing ex- 
pectation value. A general expression for the coefficients 
c a reads 

c a = \<o* V J] U iS > (10) 

r(i,j = oo) 

The correlation function has been modified by the path- 
ordered insertion of a string Y of links connecting the 
point i to 00 making it path-dependent; the modification 
ensures that the density matrix is gauge- invariant. In a 
similar fashion, we can obtain the two-particle reduced 
density matrix by expanding it in terms of the tensor 
product of Pauli matrices at the two sites under consid- 
eration: pij = j J^a /3=o CQ < acr ? ® a j- The coefficients of 
the expansion c a p are the connected correlators of two 
Wilson lines. The two-particle reduced density matrix is 
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useful in calculating the amount of entanglement localis- 
able between the two chosen spins [2||. Similar results 
hold for three-particle (as in the GHZ-state discussed by 
Arvind) and, in general, for n-particle entanglement. In 
each of these cases, the density matrix is expressed in 
terms of correlators of gauge invariant observables in a 
topological field theory i.e. topological invariants [30| . 
The difference between two choices of the path T is a 
measure of the frustrations enclosed by the loop formed 
by the two paths. There is another way in which a closed 
loop can be obtained from an open path namely, by im- 
posing periodic boundary conditions. If we choose this 
option, c a and c a p are just the expectation values of Wil- 
son loop observables in the lattice Chern-Simons theory. 

Although the Wilson loop observables are gauge- 
invariant, they are by no means the only interesting ob- 
servables. Notice that the path connecting two points on 
the primary lattice pierces one plaquette on the dual lat- 
tice with every step it advances, accumulating one unit of 
flux in the process. Let us therefore consider the gauge- 
invariant operator C = V^ 1 UV obtained by dressing 
(conjugating) the string Uij by the group-valued fields 
V on the dual lattice. Labelling the dual lattice sites 
by barred coordinates, if V runs from site i to site j, 
V^ 1 runs in the opposite direction circumnavigating the 
link U. We can use this operator instead of the Wilson 
loop to define the reduced density matrices. It may be 
mentioned that because of duality, the above conjugation 
operation simply corresponds to local unitary transfor- 
mations of nearest neighbour spins on the dual lattice. 
Physically the operator C represents a tube of dual pla- 
quettes whose axis lies on the primary lattice with fixed 
end-points. 

So far there is nothing quantum about the discussion 
of entanglement. Indeed, the Ising spins we have con- 
sidered take values ±1, much like classical bits; they are 
not allowed to be in any superposition state. The den- 
sity matrices that we considered are purely thermal in 
nature, obtained, as they are, from the Gibbs' measure. 
This kind of entanglement is called thermal entanglement 
[3l|. By using the standard Suzuki- Trotter method, 
however, we can map the three-dimensional Ising spin 
glass to a two-dimensional Ising spin glass in a trans- 
verse magnetic field. The presence of the transverse mag- 
netic field allows for transitions between the two classical 
states of the Ising spins and makes the system quantum 
mechanical. Care must be exercised in defining the spin- 
flip operation: a z and a x behave differently under gauge 
transformations because of their noncommutativity |32j . 
Furthermore, the situation here is slightly more com- 
plicated because the Suzuki- Trotter mapping from a d- 
dimensional classical statistical mechanical system to the 
(d — l)-dimensional quantum system, requires the classi- 
cal system to have different couplings along the missing 
(replica) dimension and the remaining dimensions. This 
makes the duality transformation technically a little more 



involved. The results in so far as the entanglement (now 
truly quantum in nature) are concerned follow the same 
pattern. The correlation functions that appear in this 
case are those of the quantum spins on a two-dimensional 
square sub-lattice of the original three-dimensional lat- 
tice whose third dimension acts as the discretised time 
direction. To summarise, the reduced density matrices 
that one is interested in, both for thermal and quantum 
entanglement, are expressed in terms of expectation val- 
ues of gauge invariant operators in a topological gauge 
theory. They are therefore topological invariants. 

It is difficult to obtain a ready physical insight into the 
correlation functions that appear above. Let us therefore 
consider a simpler case. Recall that similar results hold 
for all abelian groups Z p . In particular for p — > oo, we 
have a U(l) Chern-Simons gauge theory on the lattice 
and all the paths are continuous. In this case, it is well- 
known that the topological invariant produced by a com- 
putation of the correlation function of a pair of Wilson 
loops is the Gauss's linking number. The operator C we 
introduced, is the lattice generalisation of the operator 
introduced in the continuum BF theory by Cattaneo et 
al. [3^, and gives it a nice physical interpretation as a 
tube of dual plaquettes. Its correlation function gives the 
Alexander-Conway polynomial of the (possibly) knotted 
axis of a plaquette-tube, with fixed end-points on the 
primary lattice. The partition function Z is also a topo- 
logical invariant, namely, the Reidemeister torsion of the 
three manifold defined by the boundary conditions we 
choose to put on our lattice. It is equivalent to the corre- 
sponding invariant in the continuum theory, namely, the 
Ray-Singer analytic torsion 18]. Curiously enough, the 
entanglement properties of the GHZ-state which has a 
non-zero tripartite entanglement, but a zero bipartite en- 
tanglement cannot be accounted for, by the simple link- 
ing number invariants. This property is exactly like the 
corresponding property of Borromean rings which, how- 
ever, are known to be distinguished from a disjoint union 
of unlinked rings by a higher order topological invari- 
ant, namely the Massey triple product [34|. An abelian 
topological theory cannot produce this invariant. It is 
therefore necessary to treat the spins as genuinely non- 
abelian objects, like one is forced to in the presence of a 
transverse magnetic field. In view of these remarks, the 
intriguing parallels between the entanglement of quan- 
tum states and the entanglement of braids, knots and 
links, that first appeared in 0, E|, are more realistic - 
may be even natural and deep-rooted - than hitherto 
expected. 

Intuitively, topological invariants arc insensitive to the 
presence, or changes, in length scales; this being the only 
distinction between a lattice and the continuum, we could 
borrow the relevant topological invariants from the sim- 
ple continuum U{1) theory in the above discussion. Such 
a luxury is lost if the group under consideration is fi- 
nite, e.g. Z2 in the case of the Ising model. Unlike 
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classical gauge symmetries which come from continuous 
local invariances, finite groups usually appear as rem- 
nants of a continuous symmetry group which is spon- 
taneously broken. In such instances, these groups are 
known to give rise to cocycles in quantum field theory. 
This is a reflection of nontrivial group cohomology. In 
the present case, the finite group appears because we are 
dealing with spin systems. The construction of topologi- 
cal theories with finite groups follows from a deep result 
due to Dijkgraaf and Witten 0], who showed that the 
Chern-Simons actions for a finite group H are in one to 
one correspondence with the elements of the cohomol- 
ogy group H 4 (BH, Z) , BH being the classifying space 
of H. The isomorphism H 4 (BH, Z) » H 3 (H, U{1)) fur- 
ther implies that these actions are algebraic 3-cocycles, 
lu e H 3 (H,U(1)), which take values in U(l). The pres- 
ence of no-trivial group cohomology in general leads to 
non-trivial G bundles and the path integral involves a 
summation over all possible bundles. Unlike in the U(l) 
theory, in which a determinantal expression can be de- 
rived for the partition function [35j . the resulting invari- 
ants in the case of finite groups - called the Dijkgraaf- 
Witten invariants - are not expressible in terms of de- 
terminants. Analogues of linking numbers in topological 
gauge theories with finitegauge groups have also been 
worked out by Ferguson [3flj . 

As the first closing remark, we mention that the de- 
tails omitted in this letter can be found in a longer pub- 
lication under preparation. Second, we wish to point out 
that many interesting problems remain. The connections 
between quantum entanglement and more sophisticated 
topological invariants like the Jones polynomial, require 
a nonabelian generalization of the results of this letter. 
Similar investigations in two and four dimensions should 
produce interesting connections between quantum entan- 
glement and other important mathematical results in low 
dimensional topology like the intersection theory on the 
moduli space of Riemann surfaces, and Donaldson's in- 
variants respectively. Finally, it is not an exaggeration to 
say that this letter offers a mere glimpse of a new vista 
which is beginning to unfold, on the relevance of finite 
group cohomology in the studies of entanglement in spin 
glasses. We hope to dilate on these issues, in the near 
future. 

I thank M. Panero and Siddhartha Sen for some inci- 
sive comments and also A. P. Balachandran, T. R. Govin- 
darajan, I. Tsutsui, and A. Wipf for their interest and 
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